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We analyze the results of scanning near-field infrared spectroscopy performed on thin films of 
a-Si02 on Si substrate. The measured near-field signal exhibits surface-phonon resonances whose 
strength has a strong thicliness dependence in the range from 2 to 300 nm. These observations 
are compared with calculations in which the tip of the near-field infrared spectrometer is modeled 
either as a point dipole or an elongated spheroid. The latter model accounts for the antenna effect 
of the tip and gives a better agreement with the experiment. Possible applications of the near- field 
technique for depth profiling of layered nanostructures are discussed. 



PACS numbers: 68.37.Uv, 63.22.-mP 

I. INTRODUCTION 

Scattering scanning near-field optical microscopy (s- 
SNOM^pH^ is a powerful tool for probing local electro- 
magnetic response of diverse materials. The s-SNOM 
achieves spatial resolution of 10-20 nm, which is espe- 
cial ly v aluable in the physically interesting infrared re- 
gioipEl where the resolution of conventional spectroscopy 
is fundamentally limited by a rather large wavelength A ~ 
5-500 /xm. The s-SNOM techniques have been rapidly 
advancing,^ ^ which enabled their applications to imag- 
ing spectroscopy of complex oxidea^HiSl and graphene.l^ 

The s-SNOM utilizes scattering of incident light by the 
tip of an atomic force microscope (AFM) positioned next 
to the probed sample (Fig. [T]). The tip couples to the 
sample via evanescent waves of large in-plane momenta 
q ^ l/a, where a is the tip radius of curvature (a few tens 
of nm). This is why the lateral resolution of the s-SNOM 
is determined primarily by a rather than A.'^^'^ 

One of the interesting open questions is the depth {z- 
coordinate) resolution of the s-SNOM probes. Previous 
experiments suggested that it is comparable to the lat- 
eral resolution ~ a, based on imaging of small sub-surface 
particles.'^ Surprisingly, our recently near-field measure- 
ments of Si02 thin films have demonstrated that films as 
thick as several hundred nm have a response clearly dif- 
ferent from that of the bulk materiaL^Sl Thus, if instead 
of particles one has layers, then the s-SNOM is able to 
detect them at much larger depths. 

In this paper these experimental results are re-analyzed 
and compared with two theoreti cal m odels, the conven- 
tional point-dipole approximatior!-'^" and the spheroidal 
model. The former is very simple to implement but is 
also very crude. Predictably, it yields a bulk-like re- 
sponse of the s-SNOM signal as soon as the Si02 film 
thickness exceeds the tip radius, in disagreement with 
the experiment. A plausible reason for shortcomings of 
the point-dipole model is its failure to account for the 



strongly elongated shape of the tip. Such a tip acts as an 
optical antenngPt2lthat greatly enhances the electric field 
inside t he tip -sample nanogap. Unfortunately, analytical 
models,'2iJ22| ^j^g^^^ attempt to treat elongated tips do not 
apply to layered substrates. This compels us to study 
the problem numerically. 
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FIG. 1. (Color online) Schematics of an s-SNOM experiment. 
A scanned probe, modeled as a metallic spheroid with length 
2L and the apex curvature radius a, is positioned distance ztip 
above the sample. The sample contains a film of thickness di 
and dielectric function ei, which is deposited on a bulk sub- 
strate with dielectric function £2. The system is illuminated 
by infrared field f^cxt at an angle of incidence 6. Scattering of 
this radiation by the tip creates evanescent waves with large 
in-plane momenta g ~ 1 /a. The experiment measures the to- 
tal radiating dipole p of tip, which is determined by multiple 
reflections of the evanescent waves between the tip and sam- 
ple. The reflections off the sample are characterized by the 
coefficient rp{q,uj). 
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FIG. 2. (Color online) (a) Main panel: Measured infrared near-field spectra for several Si02 film thicknesses. The quantity 
plotted is the absolute value S3 of the third harmonic of the scattering amplitude normalized by that for the Si wafer, (b) 
Theoretical results for the spheroid model with a — 30 nm and L = 15a. (c) Theoretical results for the point-dipole model with 
a = 30 nm and b = 0.75a. 



To make the calculations tractable, we follow exam- 
ples in the literaturd^snssj ^^j-^j model the tip as a metallic 
spheroid of total length 2L ^ a, see Fig. [T] As shown 
below, this gives results in a much better agreement with 
the experiment in terms of both the frequency and the 
thickness dependence of the near-field signal. We at- 
tribute the origin of the more gradual film-thickness de- 
pendence in the spheroidal model to the aforementioned 
"antenna effect." The magnitude of this effect is deter- 
mined by the material response over length scales ranging 
from a to 2L, and so it truly saturates only when the film 
thickness becomes much larger than 2L. 

The remainder of the paper is organized as follows. In 
Sec. [n] we summarize the experimental procedures and 
results. In Sees. IIIII and IIVI we discuss the two theo- 
retical models and compare their predictions with the 
measurements. Concluding remarks are given in Sec. |V] 



II. EXPERIMENT 

To make the paper self-contained we summarize the 
results of our recent experiments^ in this Section. We 
investigated commercially available calibration gratings, 
which contain strips or islands of Si02 thermally grown 
on Si. The manufacturer specified thicknesses of the Si02 
layer spanned the range di = 2, 18, 22, 108, and 300 nm. 
A combination of CO2 and tunable quantum cascade 
lasers (Daylight Solutions) allowed us to cover the 
frequency range between 890 cm^^ and 1250 cm^^. The 
near-field data were collected using a Neaspec system. 

The measured s-SNOM signal represents the elec- 
tromagnetic field backscattered by the probe and the 
scanned sample. The complex amplitude s{u},t) of the 
backscattered field varies periodically with the tapping 
frequency ~ 40 kHz as the distance Ztip between the 
sample and the nearest point of the tip undergoes har- 



monic oscillations 

ztipit)^ ZQ + Azil-cosnt), (1) 

where Az — 50 nm typically. In order to suppress un- 
wanted background and isolate the part of the signal 
scattered by the probe tip, the signal is demodulated. 
Namely, we extracted the absolute values Sn(w) and 
phases 4>n{^) at tapping harmonics 

T 

.„e^^"= j |e»^*.(-,i), T^^. (2) 


The experimental results for the spectra are shown in 
Fig. [2ja). These spectra were intended to be taken in 
the tapping mode, i.e., for zero zq. However, experimen- 
tally zq can be determined only up to an additive con- 
stant ~ 1 nm. Therefore, we measured ZQ-dependence of 
S3 (the approach curves) shown in Fig. [sj^a) and selected 
the largest observed S3. Our results are in a qualitative 
agreement with previous experimental study, which re- 
ported approach curves for Si02 at a few discrete fre- 
quencies and film thicknesses. 

The data points in Fig. [2]Ja) represent the normalized 
amplitude S3(Si02)/s3(Si), where S3(Si02) and S3(Si) are 
the raw third-order demodulation signals averaged over 
the entire Si02 and Si areas, respectively. The statistical 
uncertainty of these averaged data traces is about 2%. 

For each thickness studied, the normalized amplitude 
S3(Si02)/s3(Si) exhibits several maxima. The main peak 
is situated at a; ~ 1130 cm"'^. The key aspect of the 
data is a rapid decrease in the normalized amplitude of 
this peak as the thickness is reduced. A trace of this 
resonance can be reliably identified even for the 2-nm 
thick Si02 film. Another notable feature is the growing 
strength and frequency shift of the secondary peaks on 
the high-w side of the main peak as di is decreased. 
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FIG. 3. (Color online) Approach curves, (a) Experimental data for the 105-nm thick Si02. Theoretical results for the spheroid 
(b) and the point-dipole models (c) using the same parameters as in Fig. 12] 



Since the response of Si is frequency independent in 
our experimental range, the frequency dependence of the 
spectra in Fig. [2]^ a) originates from that of Si02. We 
attribute the maxima of S3(Si02)/s3(Si) to the phonon 
modes localized at the air-Si02 interface.'^ These reso- 
nances occur in the frequency region between the bulk 
transverse and longitudinal modes of Si02 (the outer 
dashed lines in Fig. |4] below). 

The results of our theoretical calculations for the nor- 
malized scattering amplitude are presented in the re- 
maining panels of Figs. [2] and |3] They are discussed and 
compared with the experimental findings in the following 
Sections. 



III. RESPONSE FUNCTIONS AND 
COLLECTIVE MODES 

The sample is modeled as a two-layer system. The 
first layer with dielectric function ei(w) occupies the slab 
—di < 2 < 0. The second layer with dielectric function 
e2(w) occupies the half-space z < —di. The half-space 
z > { "layer 0" ) is filled with air (dielectric constant Cq = 
1). The fundamental response functions of the system are 
the reflection coefhcients rx(9,w), which are functions 
of in-plane momentum q, frequency lo, and polarization 
X = S or P. The domain of definition of rx{q,uj) is 
understood to include nonradiative modes q > y^w/c. 
It is known from previous studies that the s-SNOM signal 
is dominated by the P-polarized waves. In our two-layer 
model their reflection coefficient is given by a Fresnel-like 
formula 



where 2-axis momenta are defined by 



^2^1 — ci^i tanh ikfdi 
^ eifcf — €2^:^ tanhjfcf di 



(3) 
(4) 



Im fc| > . 



(5) 



Equation ([3| is valid for arbitrary q. In the near-field 
case where q is large and fc| ~ iq, it simplifies to 

e*-eo 62 -1- eitanhgdi 

(9, ^) - , . e* =i £1 ^ ; — ^— ^ ■ (6) 

e* 4- eo ei -|- £2 tanhgdi 



Assuming all ej are g-independent, the effective dielectric 
function {q, lu) depends on q only via the product qdi 
in this limit. Therefore, Tp {q, w) for one thickness di can 
be obtained from another by rescaling q. As discussed 
in Sec. [l] and shown in more detail below, the most im- 
portant momenta are q ^ 1/a where a '-^ 30 nm is the 
tip radius. Therefore, we can get an approximate under- 
standing of the system response by examining the behav- 
ior of Tp [q, Lj) as a function of w at fixed qdi ^ di/a. This 
behavior is dictated by the spectrum of surface collective 
modes, as follows. 

In general, surface modes correspond to poles of the 
response functions rx- Function given by E q. ([6| can 
have up to two poles at each qdi, see, e.g., Ref.1261 Thev 
are defined by the following condition on ei(a;): 



ei(w) - 



£o + £2 
2 tanh qd^ 



± 



4tanh^ qdi 



£062 ■ 



(7) 



At large qdi, where tanhgdi = 1, this condition yields 
ei{uj) = —eg or £i(a;) = —62, which correspond to modes 
localized at the upper 0-1 and the lower 1-2 interfaces, 
respectively. Actually, the latter "pole" has vanishingly 
small residue because evanescent waves do not reach the 
lower interface at qdi = 00. There is no g-dispersion 
and no coupling of the two modes in this limit. The 
dispersion appears at finite qdi, where the two modes 
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become mixed. In particular, we find 
qdi 



eo + £2 



qdi 



"0-1" 
"1-2" 



(8a) 
(8b) 



at qdi <C 1. At finite q, both, interfaces participate in 
generating these excitations. The labels "0-1" and "1- 
2" are for convenience: they indicate at which interface 
a given dispersion branch is ultimately localized as q in- 
creases. At qdi = 0, the "0-1" and "1-2" branches are 
characterized by ei(w) = and ei{uj) = —00, which cor- 
respond, respectively, to the bulk longitudinal and trans- 
verse phonon frequencies w^o and ui^^ . 

If we try to apply this formalism to real materials, we 
face the problem that Eq. ([t]) has no solutions for real 
id because the dielectric functions have finite imaginary 
parts. This is why in practice the collective mode spectra 
are usually defined differently. They are identified with 
the maxima of dissipation, i.e., Imrp. The number of 
these maxima can be fewer than the total allowed number 
of the modes because some of them can be overdamped. 
Similarly, we define w^o and as the frequencies that 
correspond to the maxima of — Ime^^(w) and Imei(cij). 

To see what kind of spectra are realized in our system, 
we use our ellipsometry data for ei(w) [Fig. [4](a)] and 
Eq. (|6| to compute Tp for several values of qdi. The 
plot of these quantities as a function of uj is presented 
m Fig.|4];c). Three maxima on each curve in the region 
of primary interest lu > 1000 cm^^ are apparent. They 
exist already at qdi — 00, and so all of them belong to 
the upper (air-Si02) interface. In fact, we do not expect 
sharp modes at the lower (Si02-Si) interface because the 
dielectric function of Si is quite large £2 ~ 11-7 in the 
studied range of w. The lowest value of Reei w —5.0 is 
not sufficient to compensate £2 and generate "1-2" modes, 
cf. Eq. 0. 

The main peak of Imrp at qdi = 00 defines the 
surface phonon frequency of Si02 uj^p ~ 1164cm^^. 
There also exist secondary peaks at cj w 1100 cm^^ and 
uj w 1220 cm~^. Their evolution as a function of qdi 
comply with the general scheme outlined above. As qdi 
decreases, all the three peaks loose strength, as expected, 
because the amount of Si02 diminishes. The lower-a; 
secondary peak redshifts, moving towards oj^^ , and then 
quickly disappears. This agrees with the Si02-Si res- 
onance being highly damped. The higher-w secondary 
peak becomes dominant at qdi < 0.5 and demonstrates 
a systematic shift towards uj^^, see Fig. 4](c). 

A notable feature of Fig. Qb) is the clustering of 
the crossing points of the different curves near lo = 
1036 cm~^. This is the frequency where the dielectric 
function of Si02 is the closest to that of Si, £2 ~ 11-7. As 
a result, the two layers act almost as one bulk material, 
so that rp(u}) is approximately thickness-independent. 

There is a qualitative correspondence between the fea- 
tures displayed by the reflection coefficient and the 
observed near- field signal S3(Si02)/s3(Si), cf. Figs. [2](a) 



and|4]jb),(c). However, the relation between rp(g,w) and 
the measured s-SNOM signal is nontrivial. For example, 
the frequency positions of the maxima in Im (g, lu) and 
those in S3(Si02)/s3(Si) differ by as much as 40cm~^. 
We also suspect that there may be some slight differences 
between the optical constants of thick films we assume in 
our calculations and those of the small Si02 structures 
we probe by the s-SNOM. This is the likely reason why 
the crossing point of the experimental curves occurs near 
1060 cm~^ rather than 1036 cm"'^ predicted by both our 
models, cf. Fig. [2j 

Developing a reliable procedure for inferring Tp (g, lo) 
from S3 remains a challenge for the theory. The next sec- 
tion presents our current approach towards this ultimate 
goal. 
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FIG. 4. (Color online) (a) Dielectric function of bulk Si02 
as a function of frequency from ellipsometry. (b) The real 
and (c) the imaginary parts of the near-field reflection co- 
efficient, rp{q,u}) for several qdi. In all the panels three 
dashed lines indicate the transverse optical phonon frequency 
'^To ~ 1074 cm~^, the surface optical phonon frequency 
ojgp « 1164 cm~^, and the longitudinal optical phonon fre- 
quency ^ 1263 cm~^. 
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IV. TIP-SAMPLE INTERACTION 

Both radiative and nonradiative waves may play sig- 
nificant roles in the s-SNOM experiment The radiative 
modes magnify the signal by a certain far-field factor 
(FFF) _F((7s,w), where — (uj/c)sm6 is the momen- 
tum of these modes for the angle of incidence 9. The 
nonradiative modes influence the effective polarizability 
x(a;,Ztip) of the tip, i.e., the ratio of its dipole moment 
and the external electric field i^oxt- Altogether the 
demodulated s-SNOM signal SnC'^'^" can be written as 



dt 

t' 



■inQ.t 



x(w,ztip(i)) 



(9) 
(10) 



Below we discuss the FFF and the tip polarizability sep- 
arately. 



Far-field factor 



The FFF for an infinite layered system is given bji^32ll 



(11) 



As shown in Figs. 5](a), for di — 300 nm Si02 film, the 
absolute value of the FFF has a maximum near cjto ~ 
1074cm^^ and a suppression near wlo ~ 1272cm~^. 
For thinner films, these features are less pronounced. The 
main maximum of S3, which is the main focus of our anal- 
ysis, is away from both a; to ^-nd wlq- H is essentially 
unaffected by the FFF. Still, if FFF were to be included 



in the calculation in the form prescribed by Eq. (Ill, it 



would produce a visible hump of 33(0;) near wto and 
a dip near wlo- These features are not present in the 
experimental data. Fig. [2]^ a). A better agreement with 
the experiment is obtained if F{qs,io) is set to a con- 
stant, which is what we do here. We rationalize this de- 
cision by noting that the Si02 layer in the actual samples 
does not extend over the entire x-y plane but occupies 
only small sub-wavelength regions. Therefore, the FFF 
is dominated by the w-independent response of Si. 



B. Point-dipole model of the tip 

The effective tip polarizability x(a;, ztip) is the most 
important factor on the right-hand side of Eq. Q and 
it is also the most difficult one to compute. This quan- 
tity is dictated by the near-field coupling between the 
tip and the sample. For irregular tip shapes it can be 
calculated only numerically. However, previous s-SNOM 
studies demonstrated that acceptable results can often be 
obtained if the tip is approximated by a spheroid ,'22H1I1 a 
small sphere j-2^'2Zi29j-i32, "finite" dipolCjEEii or a point 
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FIG. 5. (Color online) (a) The absolute value and (b) the 
phase of the far-field factor computed as a function of fre- 
quency for the incidence angle 6 = 45°. The black trace is for 
bulk Si substrate, the blue one is for bulk Si02 substrate, the 
green one is for 300 nm thick Si02 followed by bulk Si. The 
meaning of the dashed hues is the same as in Fig. |4] 



dipole.'^SlSSIlll The actual tip shape in our experiment is 
close to a rounded pyramid. 

The point-dipole approximation is the simplest one 
and it has been used extensively for modeling s-SNOM 
experim ents, including those performed on multilayer 
systems .li^lsll The point-dipole model has two adjustable 
parameters: the polarizability a"^ of the effective dipole 
and its position b with respect to the bottom of th e tip. 
The results obtained following the standard analysiJSEH 
are shown in Fig. [2](c) using a — 30 nm and b — 0.75a. 
We see that even for this rather large a the point-dipole 
model does not reproduce the observed strong depen- 
dence of S3 on thickness at di > 22 nm. 

The discrepancy can be seen more clearly in Fig. |6j 
where the height of the peak in S3(Si02)/s3(Si) corre- 
sponding to the surface phonon is plotted as a func- 
tion of c?i. For the point dipole model the curve flat- 
tens at di ^ b. In contrast, the experimentally observed 
S3(Si02)/s3(Si) maximum continues to rise with di. The 
point-dipole model also predicts a very steep approach 
curve. Fig. |3jc), in poor agreement with the measure- 
ments. 

The physical origin of the saturation of the thickness 
dependence in Fig. [2]jc) is easy to understand. One can 
think about the near-field coupling between the point 
dipole and the sample in terms of the method of images. 
For a dipole positioned at Zpd = ztip + b, the image is 
concentrated at the depth Zpd below the surface. There- 
fore, films of thickness larger than Zpd would act as a bulk 
material. Another way to arrive at the same conclusion 
is to notice that the characteristic range of momenta of 
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the relevant nonradiative waves is g ^ 1/^pd- Since 
depends on q through the term tanhgdi [Eq. ([6])], the de- 
pendence of the near-field coupling on di should saturate 
at di ^ Zpd ^ b. 



C. Spheroid model of the tip 

The lack of saturation in the observed s-SNOM signal 
as a function of c?i at di 3> a indicates that evanescent 
waves with momenta g ^ 1/a also play an important role 
in the near-field coupling between the tip and the sam- 
ple. This is a signature of models in which the tip has a 
finite extent in space 2L ^ a, see Fig. [l] Although such 
models are certainly more realistic than a point-dipole 
approximation, there has not been a systematic study of 
how the results would depend on the exact shape of the 
tip. Given some initial success of the point-dipole ap- 
proximation, we speculate that a suitable simple shape 
can provide a good compromise between increase in com- 
putational effort and ability to capture relevant physics. 

To test this idea, we model the tip as an elongated 
metallic spheroid positioned above a two-layer medium. 
This follows a tradition in the literature wherein sim- 
ilar models were considerecPsHUl fQj. (.^gg Qf bulk 
substrates. In Ref. [3T] an analytical formula for the 
spheroidal tip was also proposed, based on heuristic ar- 
guments. However, it cannot be easily extended to the 
g-dependent Tp we study here. Instead, our calculations 
are done numerically. They involve only two essential 
approximations. One is neglecting retardation, which is 
justified is the length 2L of the spheroid is smaller than 
A. The other one is neglecting the finite skin depth of 
the metal (Pt-Ir alloy) covering the tip. Due to compu- 
tational difficulties involved, this issue is left for future 
investigation. 

The calculations were performed in two ways. First is 
the standard boundary-element method. In this method 
we divide the entire tip — assuming azimuthal symmetry 
— into a large number (typically, 200) of small cylindri- 
cal segments. We assume that different segments interact 
by Coulomb interaction as coaxial rings. The interaction 
of each segment with itself is defined in such a way that 
the polarizability of the tip in the absence of the sample 
coincides with the known analytical result for the prolate 
spheroid. The effect of the sample is included by adding 
ring-ring interactions mediated by refiected electrostatic 
fields. This is accomplished by numerical quadrature over 
the product of (g, a;) and suitable form-factors. This 
is the most time-consuming step of the simulation. Af- 
ter the interaction kernel is generated in this way, it is 
straightforward to solve numerically for the dipole mo- 
ment of the tip induced by a unit external field, which is 
the desired polarizability xi^iZtip)- 

We also developed a second numerical method of com- 
puting X (to be described elsewhere) , based on an expan- 
sion of the electric field in ellipsoidal harmonics. This 
alternative method is similar to that used for a metallic 
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FIG. 6. (Color online) The thickness dependence of the S3 
peak for different tip models. The circles represent the point- 
dipole calculations, one for a = 30 nm (blue) and the other for 
a — 50 nm (green) . The diamonds are for the spheroid model 
with a = 30 nm and L = 15a, the same as in Figs.[2|^b). The 
black squares are derived from the experimental data shown 
m Fig.[2];a) after some smoothing over fluctuations. 



sphere above a dielectric half-space.'^ZlSIl We verified that 
the two methods give identical results. 

Substituting the computed polarizability x into 
Eqs. (|9| and demodulating per Eq. (10), we obtain ap- 
proach curves. Figurejsjb) illustrates that some approach 
curves are nonmonotonic near the resonances. In cal- 
culating the s-SNOM amplitude S3 we choose ztip that 
corresponds to the largest S3 because this is how it was 
done in the experiments. The results for the normalized 
amplitude are plotted in Fig.[2j^b). 

The spheroid model has two adjustable parameters: 
the apex radius of curvature a and the half-length L. 
When L = a the spheroid becomes a sphere. In this 
case the spheroid model gives results similar to the point- 
dipole model, i.e.. Fig. [2|^c). As the ratio L/a increases, 
the differences appear. However, once L/a exceeds ten, 
the normalized signal S3(Si02)/s3(Si) does not change 
much at di < 300 nm. Therefore, for long spheroids we 
effectively have only a single adjustable parameter, a. 
Remarkably, the thickness dependence of the S3 peak for 
the spheroid model matches the experiment extremely 
well (Fig.|6]). 



V. CONCLUSIONS 

In this paper we analyzed the results of experimen- 
tal study of amorphous Si02 films on Si obtained by 
scanning near-field optical spectroscopy.E^l We discussed 
the collective mode spectra of such structures and com- 
pared measurements with two theoretical calculations. 
The first is based on a conventional approximation in 
which the tip of the scanned probe is modeled as a point 
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dipole. In the second the tip is treated as an elongated 
spheroid, significantly improving agreement with the ex- 
periment. 

We explain the qualitative difference between the two 
models as follows. An important physical ingredient 
missing in the point-dipole model is the enhancement of 
the electric field near the apex of the tip — the antenna 
effect. This phenomenon is well-known from classical 
electrostatics. The enhancement of the field is controlled 
primarily by the ratio of the total length of the tip 2L 
(actually, the smaller of 2L and A) and the apex radius 
of curvature ^ a. The point-dipole model has been suc- 
cessful in the past without this enhancement factor only 
on account of the normalization procedure. Instead of 
absolute s„, one usually reports s„ normalized to some 
reference material such as Au or in our case. Si. This 
way, one eliminates any possible frequency dependence 
of the source radiation, but at the same time cancels the 
part of the signal scaling with tip size. For a stratified 
sample this cancellation is imperfect because the the field 
enhancement depends also on the dielectric response of 
the sample, which is a function of momentum q. For 
a tip of length 2L, harmonics relevant for the field en- 



hancement have momenta ranging from q ^ 1/a down 
to q ~ 1/L. Therefore, one may expect that the depen- 
dence of the s-SNOM signal on the thickness di of the 
top layer would saturate only when di ^ L. Our simu- 
lations provide direct evidence for this claim. Therefore, 
wc think that the spheroid model holds a great promise 
as an analysis tool for near-field experiments. It captures 
a lot of physics relevant to the near-field interaction while 
remaining computationally fast. 

The strong experimentally observed thickness depen- 
dence of the near-field signal^^ indicates that s-SNOM is 
capable of not only high lateral resolution but can also 
probe the system in the third dimension. However, the 
response of a layered system is different from those con- 
taining small subsurface particles^*. We hope that ex- 
perimental and theoretical approaches presented in this 
paper may be of use for accurate depth profiling of vari- 
ous dielectric and metallic nanostructures. 
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